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Abstract 

In this contribution, the evaluation of the diversity of the MIMO MMSE receiver is addressed for 
finite rates in both flat fading channels and frequency selective fading channels with cyclic prefix. It 
has been observed recently that in contrast with the other MIMO receivers, the MMSE receiver has a 
diversity depending on the aimed finite rate, and that for sufficiently low rates the MMSE receiver reaches 
the full diversity - that is, the diversity of the ML receiver This behavior has so far only been partially 
explained. The purpose of this paper is to provide complete proofs for flat fading MIMO channels, and 
to improve the partial existing results in frequency selective MIMO channels with cyclic prefix. 

Index Terms 

Diversity, Flat fading MIMO channels. Frequency selective MIMO channels. Outage probability, 
MMSE receiver 

I. Introduction 

The diversity-multiplexing trade-off (DMT) introduced by lUl studies the diversity function of the 
multiplexing gain in the high SNR regime. 121 showed that the MMSE linear receivers, widely used for 
their simplicity, exhibit a largely suboptimal DMT in flat fading MIMO channels. Nonetheless, for a 
finite data rate (i.e. when the rate does not increase with the signal to noise ratio), the MMSE receivers 
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Fig. 1. Considered MIMO system 



take several diversity values, depending on the aimed rate, as noticed earlier in 131, and also in lH, ||5l 
for frequency-selective MIMO channels. In particular they achieve full diversity for sufficiently low data 
rates, hence their great interest. This behavior was partially explained in O, iQ for flat fading MIMO 
channels and in Q for frequency-selective MIMO channels. Indeed the proof of the upper bound on the 
diversity order for the flat fading case given in |6l contains a gap, and the approach of Jll based on the 
Specht bound seems to be unsuccessfuU. As for MIMO frequency selective channels with cyclic prefix, 
fTl only derives the diversity in the particular case of a number of channel taps equal to the transmission 
data block length, and claims that this value provides an upper bound in more realistic cases, whose 
expression is however not explicitly given. In this paper we provide a rigorous proof of the diversity 
for MMSE receivers in flat fading MIMO channels for finite data rates. We also derive the diversity in 
MIMO frequency selective channels with cyclic prefix for finite data rates if the transmission data block 
length is large enough. Simulations corroborate our derived diversity in the frequency selective channels 
case. 



II. Problem statement 

We consider a MIMO system with M transmitting, N > M receiving antennas, with coding and 
ideal interleaving at the transmitter, and with a MMSE linear equalizer at the receiver, followed by a 
de-interleaver and a decoder (see Fig. [T|). We evaluate in the following sections the achieved diversity by 
studying the outage probability, that is the probability that the capacity does not support the target data 
rate, at high SNR regimes. We denote p the SNR, / the capacity and R the target data rate. We use the 
notation = for exponential equality ID, i.e. 

/(P) = / ^ lim = d, (1) 

p-^oo log p 

and the notations < and > for exponential inequalities, which are similarly defined. We note log the 
logarithm to base 2. 
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III. Flat fading MIMO channels 
In this section we consider a flat fading MIMO channel. The output of the MIMO channel is given by 



y=^/— ±lx + n, (2) 

where n ~ CN(0, Iat) is the additive white Gaussian noise and x the channel input vector, H the N x M 
channel matrix with i.i.d. entries C3sf(0, 1). 

Theorem 1: For a rate R such that log ^ 17 ^ log:;^^' with m € {1,...,M}, the outage 
probability verifies 

P(I<i?) =/j-'^(^-^+"), (3) 

that is, a diversity of 'm{N — M + m). 

Note that for a rate R < M log (i.e. m = M) full diversity MN is attained, while for a rate 
R > MlogAf the diversity corresponds to the one derived by DMT approach. This result was stated 
by 161. Nevertheless the proof of the outage lower bound in [6] omits that the event noted 23^ is not 
independent from the eigenvalues of H^H, hence questioning the validity of the given proof. We thus 
provide an alternative proof based on an approach suggested by the analysis of [2] in the case where 
R = r log p with r > 0. 

Proof: The capacity / of the MIMO MMSE considered system is given by 

M 

1 = J]log(l + /3,- 



where /3o is the SINR for the jth stream: 



[I + ^H*H 



-1 



1. 



n 



We lower bound in the first place P(/ < R) and prove in the second place that the bound is tight by 
upper bounding P(/ < R) with the same bound. 



A. Lower bound of the outage probability 

We here assume that R/M > log(M/m). In order to lower bound 
the capacity /. Using Jensen's inequality on function x i-^ log x yields 

- M 



I < M log 
= M log 



1 



M 



1+ ^H*H 



-1 



< R) we need to upper bound 

(4) 
(5) 
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We note H*H = U*AU the SVD of H*H with A = diag(Ai, . . . , Am), Ai < A2 • • • < Xm- We recall 
that the {Xk)k=i,...,M are independent from the entries of matrix U and that U is a Haar distributed 
unitary random matrix, i.e. the probability distribution of U is invariant by left (or right) multiplication 
by deterministic matrices. Using this SVD we can write 

-1 M / M ITT i2 \ -1 

M 



M 

-T 

M ^ 



^ M . M 



+ —[Xk 



(6) 



1 ) Case m = 1: In order to better understand the outage probability behavior, we first consider the 
case m = 1. In this case R/M > logM. We review the approach of lU III], which consists in upper 
bounding © by (l + ^Ai) ^ Y.f=i mV' as Efcli ^ Using this bound in © gives 

M 



I <M log 



Therefore 



1 + 



M 



_P_ 

M' 



|Ui,f 



j=i J 



In order to lower bound P(/ < i?), iH introduced the set 

1 



1 ^ 



|Uij|2 



<M + e 



for e > 0. Then, 



\I <R)>F{{I < R)nAi] 



> 



> 



2R/M 



2R/M 



where the last equality comes from the independence between eigenvectors and eigenvalues of Gaussian 
matrix H*H. It is shown in [2, Appendix A] that F{Ai) 7^ 0. Besides, as we supposed 2^/*^ > M, we 



can take e such that 
such that 



2R/M 



> 1, ensuring that 



2R/M 



M+e 



7^ 0. Hence there exists k > 



\I<R)>¥[Xi< 



K 



which is asymptotically equivalent to p' 



-{N-M+l) 



in the sense of © (see, e.g., (H Th. II.3]). 
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2) General case 1 <m < M: By the same token as for m = 1 we now consider the general case - 
we recall that we assumed that \og{M/m) < R/M. We first lower bound Ylk i+^^x which appears in 
^ by the m first terms of the sum and then use Jensen's inequality applied on x i-> x^^, yielding 



M 



u 



> 



Using this inequality in Q, we obtain that 



1 

M ^ 



1 + — H*H 

M 



-1 



Er=iiu.,iMi+i7Afc)' 



n 



,=1 (Er=iiu..f)' 



fc=i 



(7) 



where (5a;(U) = jj Xl7=i /v^i [Vt p^^ - Equation (|7]), together with ([5]), yields the following inclusion: 

ff;5,(U)(l + ^A.)<2^/*^') c(/<i?). 



, fc=l 



Similarly to the case m = 1, we introduce the set Am defined by 



M 

5a;(U) < ^ + e, k = 1,... ,m 



for e > 0. We now use this set to lower bound P(/ < R). 

P(/ <R)>F{{I < R)nAn,) 

EW) (i + ^A..)<2^/*Mnyi. 



> 



^ \ k=l 
' m 



> 



< 



, k=l 



2R/M 

'm. 



m 



.k=l 



< 



n Am 
2R/M 



M 



+ S 



The independence between eigenvectors and eigenvalues of Gaussian matrix H*H justifies the last 



equality. As we assumed that log(Af/m) < R/M, that is m < 

2^k=l "T Af^fcj ^ M/m^+e 



2-R/M 

WJm? 



, we can choose e such that 



m < 



2R/M 

TJJmF+s 



That ensures that 



7^ 0. We show in Appendix |A] that 



this probability is asymptotically equivalent to p M+m) ^j^g sense of ([T]), leading to 



P(/ <R) > 



pm{N-M+m) ' 



(8) 
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We still need to prove that F{Am) ^ 0. Any Haar distributed random unitary matrix can be 
parameterized by independent angular random variables {ai, . . . , au^) = ct whose probability 
distributions are almost surely positive (see ||9l, ifTOl and Appendix ICl). We note the functions such 
that U = <I>m(Q;). Consider a deterministic unitary matrix U^, such that |(U^,)jjp = and denote 

by cc* a corresponding Af ^ dimensional vector. It is straightforward to check that 6^ o <I>m(a*) = M/m?. 
Functions cx ^ {5^ o <^*m)(ct) are continuous at point ct* for 1 < < m and therefore there exists 
r/ > such that the ball 23 {cx^,,ri) is included in the set |a, [dk o < ^ + k = I, . . . , m}. 

We have therefore F{Am) / as 

P(^m)= / picx)da 

J {{Sk o <S>^){a)<^+e,k=l,...,m} 

> / p{a)dcx > 

J'B(a:,,rj) 

Coming back to dD, we eventually have 

< ^) > pm(N-M+m) ' 

that is the diversity of the MMSE receiver is upper bounded by m{N — M + m). 



B. Upper bound of the outage probability 

We now conclude by studying the upper bound of the outage probability, showing that m{N — M + m) 
is also a lower bound for the diversity. Note that this lower bound has been derived in IJI, ||6l using 
however rather informal arguments; we provide a more rigorous proof here for the sake of completeness. 

We now assume that R/M < log(M/(m — 1)), i.e. m — 1 < M2^^/^^ . Using Jensen inequality on 
function y i-> log(l/y), the capacity / can be lower bounded: 



I 




I + -^H*H 

M 



)"■] 









> -Mlog — Tr 



1 + — H*H 

M 



-1 



which leads to an upper bound for the outage probability: 



V <R)< 



Tr 



I + — H*H 

M 



-1 



> 



(9) 



We need to derive the probability in the right-hand side of the above inequality. Noting 23 o 
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{Ai < A2 . . . < Am, Ek=i (1 + ^Afc)"' > M2-«/A^}, 



Tr 



I+— H*H ^ 

M 



p(Ai, . . . , AM)dAi . . .(/Am- (10) 

We now introduce fim = sup(^Xi,...,Xm)&'S„{p ^^m} and prove by contradiction that < +oo. If 
fJ-m = +00, there exists a sequence (A^"\ A2"'', . . . , A^^)neN such that A^"'' — t- +oo for any k > m. 
Besides, 

fe=l k—m 

In particular M2~^/^^ < (m - 1) + I]fclm(l + l&'^i"^) ^ which, taking the limit when n ^ +oo, 
leads to m — 1 > M2^^/*^, a contradiction with the assumption m — 1 < M2~^/*^. Hence, /im < +cx3. 

We introduce the set 23i = {Ai < A2 . . . < Am, < Aa; < k = l,...,m}, which verifies 
So C Bi. Using © and ([TOll, this impUes that 



¥{I<R)< / p(Ai,...,AM)dAi...dAM, 

which is shown to be asymptotically smaller than p-^i^-M+m) jj^g sense of ([TJ in Appendix iBl The 
diversity is thus lower bounded by m{N — M + m), ending the proof. ■ 

IV. Frequency selective MIMO channels with cyclic prefix 

We consider a frequency selective MIMO channel with L independent taps. We consider a block 
transmission cyclic prefix scheme, with a block length of K. The output of the MIMO channel at time 
t is given by 

L-l 



/=o * 

where is the channel input vector at time t, ~ C3\r(0,lAr) the additive white Gaussian noise, is 
the N X M channel matrix associated to channel tap, for Z G {0, . . . , L — 1}, and H(z) denotes the 
transfer function of the discrete-time equivalent channel defined by 

L-l 

H(z) = J]H,z-'. 

1=0 

We make the common assumption that the entries of are i.i.d and CN(0, 1) distributed. We can now 
state the second diversity theorem of the paper. 
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Theorem 2: Assume that the non restrictive condition K > M'^{L — 1) holds, ensuring that log — < 



log nvr + 



(L-l)(A/-(m-l)) ' 
K 



for any m = 1 , . . . , M. Then, for a rate R verifying 



m G {1, . . . , M}, the outage probabiUty verifies 



m-l , (L-l)(Af-(m-l)) 



+ 



K 



-m(LN~M+m) 



(11) 



(12) 



P(/ < R) -- 

that is a diversity of m{LN — M + m). 

The diversity of the MMSE receiver is thus m{LN — M + m), corresponding to a flat fading MIMO 
channel with M transmit antennas and LN receive antennas. For a large block length K, the upper bound 
for rate R is close to the bound of the previous flat fading case log . Concerning data rates verifying 
- log (^5^ + ^(M - (m - 1))) < ^ < log the m{LN - M + m) diversity is only an upper 
bound; nevertheless the diversity is also lower bounded by (m — \){LN — M + {m — 1)). 



Proof: Similarly to previous section the capacity of the MIMO MMSE system is written I 



Ylif=i log(l + Pj)^ where fij is the SINR for the jth stream of x^. It is standard material that in MIMO 



frequency selective channel with cyclic prefix the SINR of the MMSE receiver is given by 

1 



/5. 



where S(z^) = Itv + {j'tl{e^''"')*U{e'^'''^). 



1, 



(13) 



A. Lower bound for the outage probability 
We assume that R/M > log(M/m). 

One can show that function A i-^ (A~^)jj, defined over the set of positive-definite matrices, is convex. 
Using Jensen's inequality then yields 



1 

K 



K 



k=l 



> 



JO 



K 



k=l 



1=0 



n 



The last equaUty follows from the fact that Y^k=i e^*''^^^'"") = Using this inequality in the SINR 
expression ( [T3] ) gives 



1 + /3, < 



L-l 



1=0 



3J 
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We now come back to the capacity / of the system; similarly to (0]), using Jensen's inequality yields 



I < M log 
< M log 



1 



M 



L-1 



j=l ^ ^ L 1=0 





-1 " 


1J 





We can now use the results of section IIII-AI by simply replacing N x M matrix H in ([S]) by LN x M 
matrix H = [H^, H^, . . . , H^^]^. They lead to the following lower bound for the outage capacity, for 
a rate R verifying R/M > log(M/m): 

P(/ <R) > ^ 



pm{LN-M+m) ' 



B. Upper bound for the outage probability 



We assume that ^ < — log ( — \ 



We first derive a lower bound for the capacity /. 



M I K \ 

^ = -E'°4^E([s{¥)]-% 

j=l \ k=\ / 

(]f^ETK[s{¥)]-')) 



> -Aflog 

The latter inequality follows once again from Jensen's inequality on function x ^ log x. 

We now analyze Tr (S(z^)^^). To that end, we write LN x M matrix H = [Hq^, . . . ,H^^]^ under 
the form H = 0(H*H)i/^ where = [0^, . . . , and 0*0 = Im- Besides, we note U*AU 

the SVD of H*H with A = diag(Ai, . . . , Aa/ ), Ai < . . . < \m- Hence, 

H(e2i^'^) = 0(e2^^'^)U*Ai/2u^ 

where &{z) = X]fco^ @iz^K Using this parametrization, 

P 



Tr (S(z^)-^) = Tr 
< Ti- 

where ^{v) = Amm(0*(e^*'''^)0(e2*'^'^)). Coming back to the outage probabiUty, 



IP(/ < R) <F 



K-l M 



fc=0 j=i ^ 



pA, 

M '\K 



> 



(14) 
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where So = {H, T.to E,^i(l + ^li^W > M2-m^}. 

We now prove by contradiction that /i^ < +00, where /x^ = supjjg^^jpAm}- If /^m = +00 there 
exists a sequence of matrices H^") gBq such that pXm^ +00. Besides, 

fc=o i=i ^ \ / / 

fc=0 j=m ^ \ / / 

As ©("^ belongs to a compact we can extract a subsequence 0('^(")) which converges towards a matrix 
0OO- For this subsequence, inequality ([TSl l becomes 

M2-f <(»-l) + iEE(l + ^7'*<"»(^)) . (16) 

k=0 j=m ^ \ / / 

Let 7oo be the function defined by 7oo(i^) = Amin(0J^(e^*'^'')0oo(e^*'^'^)) and ki, ... ,kp he the integers 
for which jooikj/K) = 0. Then det&ooiz) = det {J^fjo ®oo,iz~^) = for all z G ^e^i^''^/^,j = 
1, . . . Nevertheless, polynomial z 1— J2b=o ®oo,«-2~' has a maximum degree of M{L — 1), therefore 
p < M{L - 1). InequaUty (O then leads to 

fc^{fci,...,A:p}j=m ^ \ / / 

Moreover, ifk^{ki,.. .,kp}, X^^^'^'^^W^)) {^) ^ +00 for j > m, as 7(^("))(f ) ^ 7oo(#) 7^ for 

k ^ . . . , kp). Therefore taking the limit of ( [TT] ) when n — > +00 gives 

X Mil - 1) 
K 

which is in contradiction with the original assumption 2~^/^^ < ^^^^^ + -^^(M — (m — 1)). Hence 
/im < +00, and 23o C 23i = {H,/>Am(H*H) < prn}- Using ([HI), we thus have 

P(/ <R)< P(H G Bi), 

which, by Appendix |Bj is asymptotically smaller than p-"^i^L-M+m) jj^g sense of ([T|l, therefore 
ending the proof. ■ 

V. Numerical Results 

We here illustrate the derived diversity in the frequency selective case. In the conducted simulation we 
took a block length of = 64, a number of transmitting and receiving antennas M = N = 2, L = 2 
channel taps and an aimed data rate R = 3 bits/s/Hz. Rate R then verifies (fTTI ) with m = 1, therefore 
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SNR{dB) 



Fig. 2. Outage probability of the MMSE receiver, L=2, K=64, M=N=2 



the expected diversity is LN — M + 1 = 3. The outage probability is displayed on Fig. |2] as a function 
of SNR. We observe a slope of —10^'^ per decade, hence a diversity of 3, confirming the result stated 
in part UV] 

VI. Conclusion 

In this paper we provided rigorous proofs regarding the diversity of the MMSE receiver at fixed rate, 
in both flat fading and frequency selective MIMO channels. The higher the aimed rate the less diversity 
is achieved; in particular, for sufficiently low rates, the MMSE receiver achieves full diversity in both 
MIMO channel cases, hence its great interest. Nonetheless, in frequency selective channels, the diversity 
bounds are not tight for some specific rates; this could probably be improved. Simulations corroborated 
our results. 

Appendix A 

We prove in this appendix that, for 6 > 0, IP(E™=i P^k < b) > p-miN-M+m)_ 
We note Cm the set defined by Cm = {Ai, . . . , : < Ai < . . . < Am, J2T=i P^k < b}. As the Aj 
verify < Ai < . . . < Am, we can write 

( ™ \ r f+oo f+QO 

VpAfc<6 =/ / .../ pM,7v(Ai,...,AAf) dAi.-.dAM, (18) 

^ J J{Ai,...,A„)ee„ Ja„ Jxm-1 

where pm,n '■ I^^^ — > R is the joint probability density function of the ordered eigenvalues of a M x M 

Wishart matrix with scale matrix Im and N degrees of freedom, given by (see, e.g., |[T1): 

M 

PM,N = K-]^ n (Af -""e-^') HiX, - Xjf, (19) 

i=l i<j 
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where Km,n is a normalizing constant. We now try to separate the integral in (fTSl ) in two integrals, one 
over Ai, . . . , Am, the other over Am+i, . . . , Am- As we have (Ai, . . . , Am) £ Cm in (fTSl) . Xm < b/p and 
thus 

Pm,n{Xi, . . . , Aa/) dAm+i • • • dXM 



> / Pm,n{Xi, . . . , Xm) dXm+i ■ ■ ■ dXu 

J(A„+i,...,AM)eD 

where D = {{Xm+i, ■ ■ ■ Am) G M^"™; h/p < Xm+i < ••• < Aa/}. This integral can be 
simplified by noticing that PA/,Af(Ai, . . . , Xm) explicit expression ( fT9l ) is invariant by permutation of its 
parameters Ai, . . . , Aa/, in particular by permutation of its parameters Am+i, • • • , Aa/- Therefore, noting 
S = Sym({Am+i, • • • , Aa/}) the group of permutations over the finite set {Am+i, • • • , Aa/}, we get 



hoo f+oo 



Pm,n{Xi, . . . , Xm) dXm+i ■ ■ ■ dXM 
b/p J b/p 



/ Pm,7v(Ai, . . . , Xm) dXm+i ■ ■ ■ dXM 

sgS Js{\,„+i,...,\m)&D 

Card(S) / Pm,n{Xi, Xm) dXm+i ■ ■ ■ dXM 

J(X,„., A.,)GD 



{M-m)\ pm,n{Xi,...,Xm) dXm+i---dXA4. (21) 

'(A„+i,...,Am)GD 



Using (go]) and UJ} in ([T8), we obtain 



'( VpAfc < 6 1 > — - / 



+00 r+oo 

■ ■■ Pm,n{Xi, . . . , Xm) dXi . . . dXM- 

e„ Jb/p Jb/p 



We now replace pm,n by its explicit expression ([T9l ) and then try to separate the m first eigenvalues from 
the others. Note that we can drop the constants (M — m)! and Km,n as we only need an asymptotic 
lower bound. 

y^pX,<b]> ... n (Af -^^e-^') n(A. - A,)2 dX,... dXM 

fri J J Jb/p Jb/p ^f^) 

J&,Jb/p Jb/p Vf=l^ ' i<jlm J 

■[ n (Af-^e-^-) n (^^-^.f n {\-^,)^']dx,...dXM 

\ i=m+l i<m<j m<i<j J 
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For i < m < j, we have that Xi <b/p and thus (Aj — Aj)^ > (Aj — ^) . Hence, 

(m \ y „ m X 

^pXk<b\>( H (Af -^^e-^') n (Ai - Xjf dAi . . . dA„ J (22) 
k=l J ^"'<^'"i=l i<j<m ' 

/ • • • / n Af-^^e-^- n - ^ ) IK^^ - ^Wi . . . dAM) 



j=m+l^ ' m,<i<j 

We now have two separate integrals. We first consider the second one, in which we make the substitution 

Pi = Xi — b/ p for i = m + 1, . . . , M . 

M M / r\2m 



■■■ n (^fe-'') Ui^J--) Ui^- dXm+i . . . dXM 

Jb/p Jh/p j=m+h P' ■m<i<j 



1 p-\-oC) p-\-oo 



3m (23) 

i=m+l m<i<j 



for /9 large enough, i.e. such that e~^^^~^^^^^ > 1/2. It is straightforward to see that the integral in (l23l) 
is nonzero, finite, independent from p and therefore asymptotically equivalent to 1 in the sense of ([T||. 
Hence, we can drop the second integral in (l22l) . leading to: 



(m \ „ m 

^pXk<b\> / ll[xf-^'e-^^) H {Xi-Xjf dX,...dXm. (24) 
A;=l / i=l i<j<m 

Making the substitution ctj = pXi for i = 1, . . . , m in (l24l ) and noting Q'^ = {ai, . . . , am '■ < ai < 
■ ■ • < Om, ZlfcLi CKfc < b} we then have 

A:=l / ^ -^^'^ i=l i<j<m ' 

^ ^-miN-M+m) I -Q (^„iV-M^-«.^ _ ^^.)2 ^ ^ ^ (25) 

for /9 > 1, as we have then e~"'/^ > e~"* for i = 1, . . . , m. As 6 > it is straightforward to see that the 
integral in (|25] ) is nonzero but also finite and independent from p; it is therefore asymptotically equivalent 
to 1 in the sense of ([T]), yielding 



, k=l 



which concludes the proof. 
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Appendix B 

We prove in this section that P (Bi) < p-"i{M^N+m)^ where the set Si is defined by 
Si = {Ai < A2... < \m, < Afe < 6, A; = l,...,m}, 

with 6 > and Ai, . . . , Aa/ the ordered eigenvalues of the Wishart matrix H*H. We use the same 
approach as in Appendix |Al For we note pm,n the joint probability density function of the ordered 
eigenvalues of a M x M Wishart matrix with scale matrix and N degrees of freedom, the probability 
P(23i) can be written as 

P(Bi) = / PmA^i^ ■ ■ - ^m) d\i... d\M- 

J(Ai,...,Am)GSi 

Similarly to Appendix lAl we try to upper bound P(Si) by the product of two integrals, one containing 
the m first eigenvalues and the other the M — m remaining eigenvalues. We first replace pm,n by it 
explicit expression ([T9l ): 

. M 

P(Si) = K^l^ / H Af -^e-^' HiX^ - \,f dXi... d\M 

i(Ai,...,AM)GSi 
J(Ai,...,Am)G2Ji \^.^-^ i<j<m / 

•( n (Af-^v^-) n (^^-^.)' n o^^-^^A 

\ i=m+l i<m<j m<i<j / 



dXi . . . d\M ■ 



Note that we dropped the normalizing constant Km,n, as K^^^ = 1. For i < m < j, v/e have 
|A, - A,I < A, and thus U^<m<M^ " < U'lrn+i , yielding 

fb/p fb/p fb/p f+co r+00 / , , \ 

< /•••/ / •••/ n Af --e-' n - ^.)^ 



M \ 

n (Ar^™-^^e-^^) n (A.-A,f 

=m+l m<i<j / 



dXi . . . d\ M 



\ i=m+l m<i<j 

In order to obtain two separate integrals we discard the Am in the integral bound simply by noticing that 
Xm > 0, therefore 

/ pb/p nb/p nb/p / X \ 

P(Si) < / / . . . / n (Af-'^e-^') n (A. - A,)^ dAi . . . dX^ 

L L ■■■L n (Af -^^^-^^^e-^-) n (a. - a,)^ dx^^, . . . dx,, 



i=m+l m<i<j 
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As the second integral (in Xm+i, ^m) is nonzero, finite and independent of p it is asymptotically 

equivalent to 1 in the sense of ([U. Hence, 

rb/p fh/p fb/p "1 



JXi JX„,-i i<j<m 

We now make the substitutions ai = pXi for i = 1, . . . , m inside the remaining integral. 

rb/p fb/p fb/p rn 

/ / • • • / n n - ^.)' dx,... dXm 

Jo JXi J\^-i i^-^ i<i<m 



(26) 



-m{N-M+m) 



b rb rb ™ 



<P 



-m(N~M+rri) 



Jai Ja,n-i j^—i 
b rb rb 



i<j<m 



fO fb fO 

//•••/ n "^"^'^ n - "j)' • • • (27) 

Jo Jai Ja^^i .^^ 

as e~°*/'' < 1. The remaining integral in ( |27] ) is nonzero (6 > 0), finite and does not depend on p; 
therefore, (|27] ) is asymptotically equivalent to in the sense of ([T]). Coming back to (l26l) 

we obtain 

Appendix C 

In this appendix, we review the results of ifTOll for the reader's convenience. 
It has been shown in |9] that any n x n unitary matrix An can be written as 

1 

An-1_ 

with An-i a (n— 1) x (n— 1) unitary matrix, dn a diagonal phases matrix, that is dn = diag(e*'''i , • • • , e*'^") 
with ipi,...,ipn G [0,27r], and 0„ an orthogonal matrix (the angles matrix). Matrix 0„ can be 
written in terms of parameters 6i,...,6n € [0,^] thanks to the following decomposition: 0„ = 

Jn~l,nJn-2,n~l ■ ■ ■ Jl,2^ whcrC 



COS 9i — sin 9i 
cos 9i — sin 9i 
I„_i_i_ 

Let Ujv/ be a M X M unitary Haar distributed matrix. Then, using decomposition (1281) . 

1 



(28) 



Um = Dm(v'i)Vm(0i 



Ua/-i 



November 17, 2011 



DRAFT 



16 



with (^1 = (v9i,i,...,(^i,m) e [0,27r]A^ 6>i = . . . , ^i,a/-i) G [0, f ]*^-\ Dm(<^i) the diagonal 
matrix defined by Dm(<^i) = diag(e*'''^'^, . . . , e*"^^ *^), Vm(^i) the orthogonal matrix defined by 
Va/(0i) = Jm-i,mJm-2,m-i • • • >/i,2 and Um aM — IxM — 1 unitary matrix. Matrix Ujv/-i can 
naturally be similarly factorized. 

Similarly to |[T0|| . we can show that, in order Ua/ to be a Haar matrix it is sufficient that {(pi.i)i=i,...,M 
are i.i.d. random variables uniformly distributed over interval [0, 27r[, that 9i^i, . . . , ^i,a/-i are independent 
with densities respectively equal to (sin^i)^^"^, (sin 62)'^^'^, • • • , (sin 9m -2), 1 and independent from ^pi 
and that \Jm-i is Haar distributed and independent from ipi and 61. The proof consists in first showing, 
by a simple variable change, that if the (<^i,j)j=i,...,Af and the 6*1,1, ^i,Af-i follow the mentioned 
distributions then Da/(¥'i)Va//(0i) is uniformly distributed over the unity sphere of C^^. The proof is 
then completed by showing that if \Jm-i is a Haar matrix independent from cpi and 61 then Um is 
Haar distributed. 

Finally one can parameterize a Haar matrix Ua/ by (pi, 61 and Ua/-i. Repeating the same 
parametrization for Ua/-i we obtain that \Jm can be parameterized by the following independent 
variables 

(V'l,!) • • • I'^l.A/), (6*1,1, • • • ,^l,A/-l), (V'2,1, • • • , 9'2,A^-l), (6*2,1, • • • , 6'2,Af-2), • • • , 
{fM-2,1, ^M~2,2),()M-2,l,fM-l,l, 

whose probability laws are almost surely positive. 
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